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ABSTRACT 


LosCon,  the  software  program  developed  for  the  author’s  thesis  and  tested  at  sea, 
is  designed  to  help  the  ASW  commander  regain  tactical  control  in  a  loss  of  submarine 
contact  situation.  Persistent  detection  and  cueing  in  the  battlespace  depend  on  utilizing 
contact  reports  from  a  network  of  combatant  platform  and  offboard  sensors.  LosCon,  an 
extended  Kalman  filter-based  program  modeled  after  MTST  (Maneuvering  Target 
Statistical  Tracker),  can  integrate  the  sensor  network  very  efficiently.  Kalman  filtering  is 
a  method  of  recursively  updating  the  position  of  an  evading  target  and  accuracy  of  that 
position  using  imperfect  measurements.  Lines  of  bearing  to  the  contact  with  associated 
standard  deviation  bearing  errors  and  positions  with  their  standard  deviation  range  errors 
are  the  measurements  LosCon  uses  to  generate  an  ellipse  of  the  submarine’s  likely 
position  or  AOU  (Area  Of  Uncertainty).  LosCon  will  also  generate  an  expanded  AOU 
for  any  future  time,  allowing  commanders  to  correctly  estimate  the  size  of  the  search 
area.  The  effectiveness  of  the  sea  shield  concept  depends  on  the  ability  of  organic  forces 
to  deny  the  enemy  tactical  control  of  the  battlespace  area.  Incorporating  the  information 
generated  by  LosCon  would  assist  ASW  commanders  in  maintaining  undersea 
superiority. 
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EXECUTIVE  SUMMARY 
ASW  FUSION  ON  A  PC 


Purpose: 

LosCon,  the  software  program  developed  in  this  thesis,  is  designed  to  assist  the 
ASW  commander  regain  tactical  control  in  a  loss  of  contact  situation. 

Background: 

The  changing  strategic  environment  has  de-emphasized  traditional  cold-war 
force-on-force  tactics  in  favor  of  defeating  the  asymmetric  submarine  warfare  threat. 
Enemy  submarines  pose  a  significant  threat  to  the  Sea  Power  2 1  concept  of  operations. 
There  are  too  many  potential  adversaries  for  the  U.S.  submarine  force  to  control  the 
entire  undersea  environment,  yet  the  sea  base  and  sea  strike  forces  must  retain  their 
freedom  to  operate  without  the  interference  of  enemy  submarines.  In  the  process  of 
acquiring  undersea  supremacy,  some  periods  of  lost  contact  are  to  be  expected. 
Minimizing  these  periods  is  vital  to  maintaining  undersea  superiority. 

Discussion: 

Persistent  detection  and  cueing  in  the  battlespace  depend  on  utilizing  contact 
reports  from  a  network  of  combatant  platform  and  offboard  sensors.  LosCon,  an 
extended  Kalman  filter  based  program  modeled  after  MTST  (Maneuvering  Target 
Statistical  Tracker),  can  effectively  integrate  the  sensor  network.  Kalman  filtering  is  a 
method  of  recursively  updating  the  position  of  an  evading  target  and  accuracy  of  that 
position  using  imperfect  measurements.  Lines  of  bearing  to  the  contact  with  associated 
standard  deviation  bearing  errors  and  positions  with  their  standard  deviation  range  errors 
are  the  measurements  LosCon  uses  to  generate  an  ellipse  of  the  submarine’s  likely 
position  or  AOU  (Area  Of  Uncertainty).  The  AOU  is  updated  with  each  new  contact 
report.  LosCon  will  also  generate  an  expanded  AOU  for  any  future  time.  This  allows 
commanders  to  estimate  the  size  of  the  area  that  would  need  to  be  searched  after  an 
extended  period  of  lost  contact. 


LosCon  was  developed  following  the  observation  of  an  ASW  exercise  at  sea  with 
the  JCS  Strike  Group  in  the  Southern  California  Operating  Areas  (SOCAL).  It  was  then 
tested  at  sea  during  an  exercise  with  the  2nd  Expeditionary  Strike  Group.  Following  a 
loss  of  contact,  a  search  plan  was  designed  based  on  the  AOU  generated  by  LosCon  for 
the  time  when  maritime  patrol  craft  would  arrive  on  station.  The  aircraft  found  a 
surfaced  submarine  near  the  center  of  the  AOU. 

Conclusions: 

The  effectiveness  of  the  sea  shield  concept  depends  on  the  ability  of  organic 
forces  to  deny  the  enemy  tactical  control  of  the  battlespace  area.  Minimizing  periods  of 
lost  contact  is  vital  to  maintaining  undersea  superiority.  Incorporating  the  information 
generated  by  LosCon  would  assist  ASW  commanders  in  maintaining  undersea 
superiority.  A  further  evaluation  of  LosCon  for  operational  utilization  should  be 
considered. 
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I.  ANTI-SUBMARINE  WARFARE 


A.  SUBMARINE  THREAT 

The  twenty-first  century  has  seen  the  rise  of  a  variety  of  nontraditional  threats 
including  transnational  terrorist  organizations  and  the  rogue  nation  states  that  are  their 
traditional  sponsors.  Several  of  these  states  possess  sophisticated  asymmetric  naval 
forces,  including  submarines.  The  changing  strategic  environment  has  de-emphasized 
traditional  cold-war  force-on-force  tactics  in  favor  of  defeating  the  asymmetric  submarine 
warfare  threat.  Enemy  submarines  pose  a  significant  threat  to  the  Sea  Power  2 1  concept 
of  operations.  There  are  too  many  potential  adversaries  for  the  U.S.  submarine  force  to 
control  the  entire  undersea  environment,  yet  the  sea  base  and  sea  strike  forces  must  retain 
their  freedom  to  operate  without  the  interference  of  enemy  submarines.  In  the  process  of 
acquiring  and  maintaining  undersea  supremacy,  some  periods  of  lost  contact  are  to  be 
expected.1  “Technology  advances  and  exploitation  permit  the  pursuit  of  additional 
approaches  to  ASW  not  previously  available  [and]  sophisticated  decision  support  systems 
[...]  will  result  in  increased  probability  of  detection,”  according  to  Chief  of  Naval 
Operations  Admiral  Clark’s  ASW  CONOPS.2 3  The  software  program,  LosCon,  is 
designed  to  help  war-fighters  quickly  regain  tactical  control  in  a  loss  of  contact  situation. 

Regaining  tactical  control  is  a  complex  process.  When  loss  of  contact  occurs, 
individual  units  will  attempt  to  relocate  the  enemy  submarine,  but  the  needs  associated 
with  avoiding  counter  detection  make  the  problem  difficult.3  Tracking  a  particular 
enemy  submarine  for  an  extended  period  of  time  or  over  a  wide  area  increases  the 
difficulty  of  the  problem.  To  assist  in  the  decisions  on  how  to  conduct  the  search, 
LosCon  uses  all  the  information  available  to  generate  an  AOU  (Area  Of  Uncertainty)  for 
the  current  time  and  any  time  in  the  future  when  assets  may  become  available. 


1  Clark,  Vern.  Anti-Submarine  Warfare  Concept  of  Operations  (Draft  26  April  2004). 

2  Ibid,  7. 

3  U.S.  Navy  Publication.  NTTP  3-21.23  Submarine  Tracking  Manual.  (Washington,  D.C.:  U.S.  GPO), 
83. 
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B. 


NATURE  OF  ASW 


The  U.S.  Navy’s  current  ASW  (Anti-Submarine  Warfare)  policy  involves 
combining  information  from  a  wide  variety  of  platform-based  sensors,  fixed  land-based 
sensors,  and  intelligence  sources.  4  Anti-Submarine  Warfare  is  often  referred  to  as  an  art 
in  deference  to  the  many  educated  estimates  made  to  track  a  submarine  over  an  extended 
period. 

Effective  submarine  tracking  involves  a  detailed  understanding  of  the  local  ocean 
environment,  the  type  of  submarine  threat,  and  the  available  ASW  assets.5  Most  assets 
have  multiple  warfare  roles,  the  requirements  of  which  may  not  often  be  compatible  with 
optimal  ASW  positioning.  Extended  periods  without  a  credible  enemy  submarine  contact 
are  not  conducive  to  maintaining  an  ASW  focus  when  assets  are  also  involved  in  other 
more  active  warfare  roles.6 

Conversely,  once  there  is  evidence  of  an  enemy  submarine  in  the  area,  especially 
following  a  successful  attack  on  U.S.  forces,  there  is  a  “natural  tendency  to  confirm  any 
unexplained  sighting  as  a  submarine  and  to  assume  that  any  periscope  or  submarine 
detection  is  an  enemy  submarine.”7  This  results  in  a  large  quantity  of  false  contact 
reports  when  there  is  no  submarine  present,  and  the  possibility  of  friendly  fire  incidents 
and  collateral  damage  when  surface  or  air  forces  mistakenly  identify  ocean  life,  friendly 
or  neutral  submarines  as  enemy  submarines.  In  his  biography  Admiral  Sandy 
Woodward,  the  United  Kingdom’s  battle  group  commander  during  the  Falklands  War, 
commented  on  frequent  “submarine  contacts”  that  turned  out  to  actually  be  whales.8  In 
WWII,  U.S.  surface  forces  sank  two  U.S.  submarines  after  mistaking  them  for  enemy 
combatants.  The  number  of  false  contact  reports  associated  with  submarine  tracking 
remains  high,  and  must  be  considered  in  any  search  algorithm  designed  to  locate  enemy 
submarines. 


4  U.S.  Navy  Publication.  NWP  3-21  Navy  ASW.  (Washington,  D.C.:  U.S.  GPO),  29. 

5  Ibid. 

6  Ibid,  34-35. 

7  Ibid,  35-36. 

8  Woodward,  John.  One  Hundred  Days:  the  Memoirs  of  the  Falklands  Battle  Group  Commander. 
(Annapolis:  Naval  Institute  Press,  1992),  96-91 . 
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A  search  for  a  submarine  often  begins  with  theater  and  or  national  intelligence. 
Local  oceanographic  infonnation  must  then  be  analyzed  to  detennine  the  best  way  to 
prosecute  the  search  considering  the  assets  available.9 

ASW  assets  encompass  a  broad  range  of  systems  that  can  provide  a 

serious  threat  to  an  enemy  submarine  while  it  is  submerged.  ASW  capable 

U.S.  assets  include: 

•  National  sensor  systems 

•  IUSS  (Integrated  Underwater  Surveillance  System) 

•  ASW  capable  surface  ships 

•  ASW  fixed  wing  aircraft  (P-3,  S-3) 

•  ASW  helicopters  (SH-60) 

•  Submarines.10 

All  of  these  assets  have  non-ASW  roles  as  well.  Thus  it  is  not  surprising  that, 
whether  attempting  to  prosecute  a  real-world  or  exercise  problem,  ASW  forces 
experience  greatest  difficulty  maintaining  tactical  control  when  the  enemy  submarine  has 
crossed  into  a  new  area  of  responsibility.11  It  may  not  be  possible  for  the  new 
commander  to  assign  sufficient  assets  to  tracking  the  enemy  submarine. 


C.  OVERVIEW  OF  THESIS 

The  effectiveness  of  the  sea  shield  concept  depends  on  the  ability  of  organic 
forces  to  deny  enemies  tactical  control  in  the  areas  occupied  by  sea  strike  groups  and  sea 
bases.  Control  of  the  sea  in  the  Mahanian  sense  is  not  feasible  without  sufficient  U.S. 
submarines  to  match  potential  enemy  submarines  one  for  one.  Loss  of  contact  on  a 
hostile  or  potentially  hostile  submarine  threatens  the  US  force’s  ability  to  successfully 
conduct  their  missions. 

9  U.S.  Navy  Publication.  NWP  3-21  Navy  ASW.  (Washington,  D.C.:  U.S.  GPO),  30. 

10  Ibid,  34-35. 

11  Ibid,  29-30. 
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LosCon  is  designed  to  assist  in  regaining  tactical  control  in  a  loss  of  contact 
situation.  The  ellipse  of  the  submarine’s  probable  location,  its  AOU,  can  be  used  to 
move  vulnerable  assets  away  from  the  threat  and  to  focus  searchers  efforts  to  most 
effectively  regain  contact. 

LosCon’ s  development  included  observing  an  ASW  exercise  conducted  by  the 
JCS  Carrier  Strike  Group  and  testing  of  the  program  in  an  ASW  exercise  with  the  2nd 
Expeditionary  Strike  Group.  LosCon’s  mathematical  core  is  an  iterated  extended  Kalman 
filter  employing  the  MTST  (Maneuvering  Target  Statistical  Tracker).  12  The  following 
chapters  will  explain  iterated  extended  Kalman  filters,  their  application  in  LosCon,  and 
LosCon’s  attributes,  and  potential  ways  in  which  LosCon  might  be  expanded. 


12  Wagner,  Daniel  H.,  “Naval  Tactical  Decision  Aids,”  Military  Operations  Research  Lecture  Notes, 
(September  1989). 
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II.  THEORECTICAL  BASIS  OF  LOSCON 


A.  ITERATED  EXTENDED  KALMAN  FILTERING  THEORY 

Kalman  filtering,  an  algorithm  frequently  used  in  localization  and  tracking  aids,  is 
a  method  of  recursively  updating  the  mean  of  a  system  based  on  a  series  of 
measurements.  The  system  is  modeled  in  two  parts:  measurement  and  movement. 

In  LosCon,  the  target’s  true  state  X  consists  of  two  location  components  and  two 
velocity  components  as  shown. 


X  = 


Cartesian  Coordinate  from  Longitude 
Cartesian  Coordinate  from  Latitude 
Velocity  in  the  East-West  Direction 
Velocity  in  the  North-South  Direction 


The  true  state  cannot  be  known  exactly,  so  it  is  treated  as  a  random  variable. 
Specifically,  X  is  assumed  to  be  a  multivariate  normal  random  variable  with  mean  p  and 
covariance  matrix  E,  symbolically  X~N(p,X).  The  symbol  ~  will  frequently  be  used  to 
indicate  the  probability  distribution  shown  on  the  right-hand  side.  The  Kalman  filter 
repeatedly  revises  p  and  X  to  account  for  the  passage  of  time  (X  changes  with  time),  as 
well  as  for  measurements. 


The  new  state  of  the  system  after  the  passage  of  time,  X ,  is  modeled  by  the 
product  of  the  movement  matrix  <f>  and  the  old  state  X  ,  summed  with  the  error 
distribution  associated  with  movement,  W  :  N(juw,Q)  . 13 

Movement  Model: 

X'  =  <f>X  +  W 

It  follows  that  X' :  N(/u',Y!)  ,  where  /u’ =  0/u  + /uw  and  Z'  =  ^Z^r  +  C-  If  ^ 
depends  onJ,  the  Kalman  filter  is  extended. 


13  Washburn,  A.  “A  Short  Introduction  to  Kalman  Filters,”  (NPS  2004),  1-4. 
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The  measurement  Z  is  the  product  of  a  measurement  matrix  H  and  the  state  of  the 
system  X,  summed  with  the  error  distribution  associated  with  that  type  of  measurement, 
V  :  N(!uv,R)M 

Measurement  Model: 

Z  =  HX  +  V 

When  Z  is  given,  it  follows  that  X:  N(/i,£),  where 
ju  =  ju  +  K(Z-juv-Hju)  and  ±  =  (/  -  KH)  I .  Here  K  =  1HT(HZ,HT  +R)~'  is 
the  Kalman  gain.  If//  depends  on  X,  as  it  does  in  LosCon,  the  Kalman 
fdter  is  extended. 

Kalman  gain  is  the  amount  by  which  the  old  estimate  for  the  state  of  the  system  is 
changed  by  introducing  a  new  measurement.  The  shock,  S  =  Z-/uv-Hju,  is  the 

difference  between  the  measurement  and  what  the  Kalman  filter  expected  the 
measurement  to  be  based  on  all  past  measurements.  Recall  that 

ju  =  ju  +  K[Z-juv-H/j);  ju  is  updated  by  summing  the  product  of  the  Kalman  gain  and 
the  shock  with  the  previous  // . 1 5  Dimensionless  shock  for  a  specific  measurement  i, 
DSi  =  S  '  ( //  X,  ///  +  R  j  St,  can  be  thought  of  as  normalized  shock,  and  it  is  used  to  tell 
when  the  shock,  St ,  is  too  large.  16 

Recall  that  H  is  the  measurement  matrix  that  defines  the  measurement’s  relation 
to  the  state,  and  <f>  is  the  movement  matrix  that  defines  how  the  state  changes  with  time. 
If  either  H  or  ^  depends  on  x,  the  respective  measurement  or  movement  model 
becomes  nonlinear.  Kalman  filtering  is  a  linear  process,  so  nonlinearity  introduces  a 
complication.  When  either  of  the  models  must  be  linearized,  the  process  is  referred  to  as 
an  extended  Kalman  filter. 


14  ibid. 

15  Balakrishnan,  A.V.,  Kalman  Filtering  Theory’,  (New  York:  Springer- Verlag,  1984),  76-96. 

16  Washburn,  A.  “A  Short  Introduction  to  Kalman  Filters,”  (NPS  2004). 
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Iteration  is  useful  in  extended  Kalman  filters  for  minimizing  error.  1 7  One  such 
case  is  when  measurements  can  be  bearings-only  as  they  are  in  LosCon.  For  a  bearings- 
only  measurement,  H  depends  onl.  Greater  detail  for  this  case  will  be  provided  in 
section  “B:  Implementation  in  Excel.”  Recall  that  in  the  measurement  model, 
jj.  =  ju  +  K(Z  -  nv  -  Hju) .  Some  value  must  be  entered  for  the  initial  estimate  of  the  mean  at 
the  time  of  the  first  measurement.  The  value  of  the  state  of  the  system  at  the  time  of  the 
first  measurement  is  referred  to  as  X2=1 :  ,V( )  where  the  index  z  refers  to  the  time 
associated  with  the  Z'h  measurement.  Let  i  be  the  index  of  the  iteration,  which  is 
distinct  from  the  measurement’s  time  index  z .  The  mean  of  the  state,  ju ,  can  be 
recalculated  repeatedly  as  shown:  Jui  =  Jui_l+  K(Z - nv  - until  » //, .  The 
iteration  converges  when  ~  //, .  Let  n  refer  to  the  number  of  iterations  at  which 
convergence  occurs.  Iteration  is  done  to  minimize  the  error  introduced  by  the  initial 
estimate  of  the  mean,  /./2=u=l .  When  iteration  results  in  convergence,  //rlM  is  replaced  by 

Mz=\,  i=n  • 

The  mean  of  the  state  does  not  always  converge,  because  in  some  cases  the 
linearized  measurement  matrix,  H  ,  will  cause  the  Kalman  gain,  K  =  Y.HT {HY.H7  +  R)~' , 
to  update  //  in  the  wrong  direction.  Divergence  during  the  iteration  process  can  occur 
when  the  initial  estimate  of  X,  //.=1  w ,  is  poor  or  the  measurements,  Z  ,  have  large  errors  . 

The  iterated  extended  Kalman  filter  is  a  useful  tool  in  spite  of  the  need  to  linearize 
non-linear  models,  because  iteration  can  minimize  the  error  introduced  in  the 
linearization  process.  The  potential  for  divergence  during  iteration  is  tolerable  because 
the  filtering  process  can  be  restarted  easily  by  selecting  a  new  initial  estimate. 
Divergence  also  serves  as  a  warning  for  the  existence  of  bad  data. 


Bell,  Bradley  M.,  and  Fredrick  W.  Cathey.  “The  Iterated  Kalman  Filter  Update  as  a  Gauss-Newton 
Method,”  IEEE  Transactions  on  Automatic  Control  38,  no.  2  (February  1993):  294-297. 
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B.  IMPLEMENTATION  IN  EXCEL 

LosCon  is  modeled  after  MTST,  the  program  with  which  Daniel  H.  Wagner  and 
Associates  pioneered  using  Kalman  filters  to  track  submarines  in  1980. 18  An  explanation 
of  LosCon’s  iterated  extended  Kalman  filter  will  follow.  When  naval  terminology  with 
similar  meanings  to  the  precise  mathematical  terminology  exists,  it  will  be  included. 

Recall  that 


Estimated  Mean  Cartesian  Coordinate  from  Longitude 

A Ion 

Estimated  Mean  Cartesian  Coordinate  from  Latitude 

Ato 

Estimated  Mean  Velocity  in  the  East-West  Direction 

Ae-W  Vel 

Estimated  Mean  Velocity  in  the  North-South  Direction 

_  An-S  Vel  _ 

The  Kalman  filter  begins  with  an  estimate  of  nz=l  i=l  and  a  covariance  matrix,  Zz=1 , 

estimating  the  accuracy  of  that  mean.19  LosCon  uses  an  extended  Kalman  filter,  because 
some  of  the  measurements  are  bearings  only.  Recall  that  for  extended  Kalman  filters 
iterating  //z=u=l  can  minimize  error. 20  The  initial  covariance  matrix  is  not  iterated  in 

LosCon  because  the  movement  model  is  linear.  Since  Zz=1  is  not  affected  by  the  iteration 


of  // ,  Zz=1  will  be  written  as  X  •  The  initial  covariance  matrix  has  the  default  value 
shown  below. 


2,= 


100,000 

0 

0 

0 


0 

100,000 

0 

0 


0 

0 

100 

0 


0 

0 

0 

100 


See  Appendix  B:  Visual  Basic  Code  to  modify  the  initialization  of  X! . 

LosCon  generates  AOUs  based  on  the  position  components  of  //  and  the  standard 
deviations  from  X  .21 

18  Wagner,  Daniel  H.,  “Naval  Tactical  Decision  Aids,”  Military  Operations  Research  Lecture  Notes, 
(September  1989). 

19  Balakrishnan,  A.V.,  Kalman  Filtering  Theory,  (New  York:  Springer- Verlag,  1984),  76-96. 

20  Kerr,  Thomas  H..  “Streamlining  Measurement  Iteration  for  EKF  Target  Tracking,”  IEEE 
Transactions  on  Aerospace  and  Electronic  Systems  27,  no.  2  (March  1991):  408-421. 

21  Washburn,  A.  “A  Short  Introduction  to  Kalman  Filters,”  (NPS  2004),  4. 
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A  measurement  in  LosCon  corresponds  to  a  contact  on  a  submarine.  Recall  that 
the  measurement  matrix,  H ,  is  a  description  of  what  kind  of  information  the 
measurements  contain.  The  value  of  H  depends  on  the  type  of  measurement. 
Measurements  in  LosCon  are  of  two  types. 

1)  Position  measurements  consist  of  the  latitude  and  longitude  of  the  target,  a 
single  standard  deviation  error  in  nautical  miles,  and  the  time  of  the 
contact.  For  these  measurements,  H  is  a  fixed  value. 


H  = 


1 

0 


0  0  0 
1  0  0 


2)  Bearing  measurements  consist  of  a  bearing  taker’s  latitude  and  longitude, 
a  bearing  on  a  target,  the  single  standard  deviation  bearing  error,  and  the 
time  of  the  contact.  H  varies  based  on  the  state. 


H  = 


Mla,~Z, 


lat 


( l1  Ion  Z Ion  ) 


( I-1  lat  Z lat )  +(/hon  Z Ion)  {M lat  Z lat)  +  (/Cn  Z Ion) 


0  0 


where  julat  is  the  “y”  component  of  the  initial  estimate  which  corresponds 
to  latitude,  /Jlon  is  the  “x”  component  of  the  initial  estimate  which 
corresponds  to  longitude,  Zlat  is  the  bearing  taker’s  latitude  converted  into 
Cartesian  coordinates,  and  Zlon  is  the  bearing  taker’s  longitude  converted 

into  Cartesian  coordinates.  Figure  1  shows  the  relationship  between  the 
four  terms. 
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Mlon  Z Ion 


\ 


Z Ion 

Mlon 

cos  6 

Mint  —  ^lat 

r  (Mia, 

^ lat  )  (  f^lon  ^ Ion  ) 

-sin  6 

( Ion  ~^lon) 

''  ( Mlat  Z ia, )  +(/Jlon  % Ion) 


Figure  1:  The  terms  in  the  measurement  matrix  of  a  bearing 

measurement  are  shown  in  trigonometric  form. 


Since  H  changes  based  on  the  /ulat  and  ju/on  components  of  the  state  for  bearing 

measurements  and  //  is  a  fbnction  of  //z=l  =  [//,„„  fi,a,  mE-w velocity 
recursively  calculating  //7=]  can  minimize  the  error  introduced  in  the  extended  Kalman 
filter  by  the  initial  estimate.22  Recall  that  in  the  course  of  this  iteration  it  is  also  possible 
for  the  values  of  //  to  diverge.  This  occurs  when  the  extended  Kalman  filter  updates  // 
in  the  wrong  direction. 


22  Kerr,  Thomas  H..  “Streamlining  Measurement  Iteration  for  EKF  Target  Tracking,”  IEEE 
Transactions  on  Aerospace  and  Electronic  Systems  27,  no.  2  (March  1991):  408-421. 
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After  measurements  are  entered,  the  filter  should  be  iterated.  If  //  converges 
within  the  first  several  iterations,  an  AOU  is  returned,  and  LosCon  has  performed  as 
desired.  If  that  does  not  happen,  shock  can  be  used  to  help  diagnose  what  went  wrong. 
Recall  that  the  shock,  S  =  Z  -/iv-Hju,  is  the  difference  between  the  measurement  and 
what  the  Kalman  filter  expected  the  measurement  to  be,  based  on  all  past  measurements. 
Dimensionless  shock,  DSj  =  S'  ( //  X,  HT.  +  R  j  Sf ,  is  used  to  tell  when  the  shock,  S', ,  is 

too  large.  23  High  dimensionless  shocks  will  occur  when  measurements  with  large  errors 
are  entered  in  LosCon  ,  or  when  the  initial  estimate  for  the  location  and  velocity  of  a 
target  submarine  is  poor.  The  challenge  of  dealing  with  high  dimensionless  shocks  is  that 
while  they  may  indicate  that  the  new  contact  report  is  false,  they  may  instead  mean  that 
the  last  few  reports  were  false  and  that  the  current  contact  report  is  true. 

ASW  remains  an  art.  LosCon  is  a  tool  to  help  the  war-fighter  detennine  which 
contact  reports  are  true  and  where  the  submarine  would  be  if  contact  reports  in  a 
particular  group  were  assumed  to  be  true.  The  toggle  switch  in  LosCon  allows  the  war¬ 
fighter  to  quickly  check  multiple  groups  of  contact  reports  against  each  other  to  find 
which  sets  could  be  accurate. 


23  Washburn,  A.  “A  Short  Introduction  to  Kalman  Filters,”  (NPS  2004). 
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III.  TACTICAL  UTILITY  OF  LOSCON 


A.  OBJECTIVE 

LosCon  utilizes  the  well-understood  mathematical  principles  of  Kalman  filtering 
to  make  tracking  submarines  in  a  localized  area  easier.  Kalman  filtering  is  commonly 
used  for  problems  that  include  significant  error  factors.  In  the  submarine-tracking 
problem,  the  error  types  are  two-fold: 

(1)  A  large  proportion  of  the  contact  reports  can  be  assumed  to  be 
false,  although  exactly  which  ones  are  which  is  not  immediately 
clear.  Sometimes  later  observations  reveal  some  reports  to  be 
impossible,  or  at  the  conclusion  of  an  exercise  the  submarine’s 
track  can  be  compared  to  the  contact  reports. 

(2)  Each  true  contact  report  also  includes  a  certain,  non-constant 
amount  of  error.  A  sonar  contact  report  will  have  an  associated 
bearing  error.  A  periscope  sighting  will  have  bearing  and  range 
error.  All  other  contact  reports  such  as  those  associated  with 
sonobuoy  fields  or  IUSS  will  have  their  own  errors  as  well. 

LosCon  minimizes  the  effects  of  false  contact  reports  and  can  use  the  error  associated 
with  true  contact  reports  to  create  an  AOU  for  the  location  of  the  enemy  submarine. 

The  most  effective  method  for  false  contact  reports  is  to  ignore  them.  Since  the 
operator  cannot  know  for  certain  which  contact  reports  are  false,  LosCon  includes  a 
toggle  switch  that  allows  the  operator  to  selectively  ignore  and  include  contacts,  running 
the  program  multiple  ways  to  visually  understand  each  contact  reports’  effect  on  the 
AOU. 

Operators  are  faced  with  a  complex  analytical  problem  after  losing  contact.  They 
have  a  number  of  old  bearings  on  the  contact  which  are  assumed  to  be  correct,  certainly 
contain  some  error,  and  may  in  fact  not  even  be  all  true.  Based  on  the  contact  reports  that 
the  operator  has  selected  as  true,  an  AOU  is  generated  by  LosCon.  As  new  information 

about  the  location  of  the  submarine  is  entered,  the  AOU  changes  and  the  operator  can 

13 


maintain  a  clear  tactical  picture  by  considering  all  the  information  in  one  view  screen.  In 
the  case  where  assets  need  to  fulfill  a  non-ASW  role  for  a  certain  time  period,  the 
operator  can  enter  that  future  time  to  see  the  expanded  AOU  for  that  future  time  based  the 
assumption  that  no  new  contact  reports  will  be  made  in  the  intervening  hours. 

LosCon  can  track  up  to  three  submarines.  Initial  contact  may  not  include 
sufficient  information  to  correctly  distinguish  the  enemy  submarines,  but  the  input  data 
remains  on  screen  and  accessible  throughout  the  exercise  so  that  an  operator  can 
redesignate  bearing  and  position  reports  as  contacts  on  a  different  submarine  or  false 
contacts  at  any  time. 


B.  DEVELOPMENT  AND  TESTING 

LosCon  was  designed  for  use  in  a  strike  group’s  ASW  command  center  to  track 
and  reacquire  submarine  contacts.  During  the  early  development  stage,  the  author  went 
to  sea  on  the  USS  JOHN  C.  STENNIS  (CVN  74),  a  Nimitz  class  aircraft  carrier,  to 
observe  an  ASW  battle  problem  that  the  carrier  strike  group  was  conducting  as  part  of  her 
pre-deployment  exercises.  While  the  carrier  has  too  much  ship  noise  to  be  an  effective 
anti-submarine  warfare  platfonn,  the  ASW  commander  is  on  board  the  carrier,  and  so  all 
the  contact  reports  are  sent  back  to  the  carrier. 

The  ASW  module  on  the  carrier  proved  to  be  crowded  with  computers  and 
people.  The  ASWO  (ASW  Officer)  was  equipped  with  computer  systems  predicting  the 
group’s  sonar  coverage,  and  also  with  chat  connections  to  many  of  the  other  platforms  in 
the  group.  Once  contact  was  made,  it  was  plotted  on  a  traditional  chart  and  the  ASWO 
was  left  to  estimate  where  the  submarine  might  be  next.  When  no  details  about  the 
accuracy  of  the  contact  report  position  were  mentioned,  sometimes  a  quarter  was 
centered  on  the  plot  and  used  to  draw  a  circle  around  the  point.  Without  the  infonnation 
needed  to  make  a  more  educated  estimate  of  the  error  there  was  no  reason  to  do 
otherwise. 

The  exercise  led  to  several  observations  about  the  on  board  ASW  operational 
analyses  process  in  use.  The  ASWO  was  not  getting  all  the  infonnation  available. 
Physical  space  in  the  ASW  module  was  limited,  and  there  was  not  a  computer  that  could 
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be  solely  devoted  to  a  new  program.  Most  importantly,  when  loss  of  contact  occurred, 
the  ASWO  did  not  have  tools  to  assist  in  refining  the  situation  for  the  purpose  of 
regaining  contact. 

As  such,  during  the  next  several  months  of  development,  LosCon’s  design  was 
refined  to  optimize  its  practicality.  LosCon  was  subsequently  tested  onboard  USS 
MOBILE  BAY.  The  Expeditionary  Strike  Group  had  an  ASW  battle  problem  during  its 
final  battle  problem  in  which  a  U.S.  nuclear  submarine  played  the  role  of  a  diesel 
submarine.  During  the  ASW  component  of  the  two-day  battle  problem,  the  search 
platforms  made  contact  with  the  submarine  and  maintained  contact  for  a  few  hours  early 
in  the  problem.  The  searchers  lost  contact  for  several  hours  but  had  maritime  patrol 
aircraft  coming  on  station  with  sonobouy  laying  capabilities.  Based  on  LosCon’s 
predicted  AOU  for  the  time  the  aircraft  would  arrive,  a  search  pattern  was  laid  out.  Less 
than  a  half  hour  into  the  search,  the  aircraft  located  the  submarine  near  the  center  of  the 
predicted  AOU. 


C.  AOU  GENERATION 

Based  on  past  contact  reports,  LosCon  seeks  to  generate  AOUs  containing  a 
submarine’s  most  likely  position  at  a  given  time.  A  series  of  figures  will  be  used  to  show 
the  connection  between  the  standard  deviation  errors  associated  with  each  contact  report 
type  and  the  generated  AOUs. 

Figure  2  shows  hypothetical  concurrent  contact  reports  to  demonstrate  the  use  of 
two  bearing  contact  reports  with  their  associated  standard  deviation  bearing  errors  to 
generate  a  set  of  current  AOUs.  Elongated  along  the  direction  of  the  lines  of  bearing,  the 
AOUs  shown  include  the  target  submarine  with  probabilities  of  39%,  86.5%,  and  99%, 
which  correspond  to  one,  two,  and  three  standard  deviations  respectively.  The  standard, 
which  LosCon  uses,  is  the  86.5%  probability  AOU. 
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Figure  2:  Lines  of  bearing  with  associated  errors  and  the  resulting 

AOUs. 


If  instead  of  lines  of  bearing,  the  contact  reports  on  the  target  submarine  were 
position  reports  from  another  unit,  the  AOUs  would  be  circular  as  shown  in  Figure  3. 


Figure  3:  AOUs  generated  using  position  reports  with  circularly 

symmetric  errors. 


Generally,  by  the  time  a  submarine  has  been  tracked  for  a  while  and  then  lost, 
there  are  multiple  types  of  contact  reports  for  the  target,  some  of  which  are  considered 
false,  and  the  generation  of  an  AOU  becomes  significantly  more  complicated.  Figure  4 
shows  such  a  scenario. 
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Figure  4:  Multiple  lines  of  bearing  and  positions  with  associated  errors 

are  combined  to  create  39%,  86.5%,  and  99%  AOUs. 

Once  a  contact  is  lost,  LosCon  can  continue  to  update  the  AOUs  for  the  current 
time  or  even  some  future  time.  Figure  5  shows  the  86.5%  AOU  at  the  time  of  contact 
loss.  In  Figure  6,  LosCon  uses  all  the  information  on  the  contact’s  position  used  in 
Figure  5  to  create  the  AOU  where  search  efforts  can  expect  to  regain  contact  for  a  future 
time. 


Figure  5:  A  standard  (86.5%)  AOU  at  the  time  of  last  contact. 
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Figure  6:  The  same  (86.5%)  AOU  is  enlarged  two  hours  later  with  no 

further  contacts. 

The  method  by  which  LosCon  converts  contact  reports  into  AOUs,  such  as  shown  in  the 
last  five  figures,  is  an  iterated  extended  Kalman  filter. 


D.  LIMITATIONS  AND  STRENGTHS 

There  are  limitations  inherent  in  using  LosCon’ s  Kalman  filtering  search  process. 
The  most  significant  limitation  is  the  need  for  an  initial  estimate  for  where  the  enemy 
submarine  could  be.  As  explained  in  the  previous  chapter,  this  cannot  be  avoided  when 
using  Kalman  filters.  A  bad  initial  estimate  for  the  enemy  submarine  starting  point  at 
time  zero  can  result  in  LosCon  failing  to  return  a  usable  result. 

A  lesser  limitation  is  the  need  for  an  estimate  of  the  target  submarine’s  long-term 
average  velocity,  /uw .  This  is  the  average  of  all  courses  and  speeds  a  submarine  could  be 

reasonably  expected  to  follow.  When  all  courses  and  all  potential  speeds  are  equally 
likely,  the  long-term  average  velocity  is  zero.  Absent  any  other  information,  zero  is  the 
logical  estimate.  This  is  the  default  setting  in  LosCon.  An  example  for  which  this  would 
be  a  correct  assumption  is  a  submarine  attempting  to  target  a  sea  base.  The  submarine 
would  likely  move  around  the  area  at  varying  courses  and  speeds,  but  no  particular 
direction  is  more  favored  than  another.  A  submarine  traveling  from  one  port  to  another 
should  not  have  a  zero  estimate.  While  it  will  vary  course  and  speed,  it  is  going 
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somewhere  and  thus  must  have  an  overall  non-zero  velocity.  If  no  credible  intelligence  is 
available  concerning  a  submarine’s  intentions,  a  zero  estimate  for  long-term  average 
velocity  is  best.  Zero  is  the  average  of  all  possibilities,  and  the  effects  of  non-zero 
velocities  on  the  center  point  of  the  AOU  are  small  compared  to  the  size  of  the  AOU.  To 
change  LosCon  to  allow  the  long-tenn  average  velocity  to  be  varied  from  zero,  refer  to 
the  comments  in  Appendix  B:  Visual  Basic  Code. 

The  openness  of  LosCon’ s  database  makes  these  limitations  tolerable.  At  any 
time  the  operator  can  select  a  new  initial  estimate  for  the  location  and  velocity  of  a 
particular  target  submarine  without  having  to  reenter  anything  else.  Similarly,  each 
contact  report  can  be  toggled  false  or  true  at  any  point,  and  each  contact  report  can  be 
marked  as  a  contact  report  on  submarine  one,  two  or  three.  All  the  report  data  entered 
remains  clearly  visible  to  the  operator  on  a  master  spreadsheet.  The  initial  estimates  for 
the  location  and  velocity  of  each  submarine  at  time  zero  along  with  the  future  time  at 
which  an  AOU  is  desired  are  located  on  supporting  spreadsheets  for  each  potential  enemy 
submarine. 

A  number  of  LosCon’s  attributes  were  designed  specifically  for  ease  of  use.  For 
example  the  program  file,  LosCon.xls,  is  small  enough  to  fit  on  an  1 .44MB  floppy  disk. 
Microsoft  Excel,  software  already  installed  on  most  ships’  computers,  was  chosen  as  the 
interface.  The  nature  of  the  Kalman  filter  algorithm  makes  the  program  operation  low 
memory  so  the  system  requirements  are  not  an  issue.  LosCon  can  be  left  as  an  open 
window  in  the  background  as  the  computer  was  used  to  do  other  operations.  The 
database  remains  open  to  changes.  Based  on  the  limitations  of  Excel  and  a  desire  to  keep 
LosCon  floppy  disk  sized,  the  number  of  submarines  was  limited  to  three.  To  account  for 
the  localization  math  being  done  before  the  contact  infonnation  was  sent  to  the  ASW 
command,  an  option  of  entering  the  contact  reports  as  positions  with  circularly  symmetric 
standard  deviations  in  nautical  miles  was  included.  The  expected  contact  report  fonn  of  a 
platform  location,  bearing  to  the  target,  and  bearing  error  was  kept  as  an  option. 

In  summary,  the  Kalman  filter  algorithm  recursively  calculates  an  AOU  for  the 
target  submarine  using  repeated  measurements  of  the  submarine’s  location.  An  initial 
state  must  be  estimated,  but  the  associated  covariance  matrix  is  set  by  default  to  be  very 
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large,  and  that  estimate  does  not  affect  the  AOU  appreciably  after  two  or  more  contact 
reports  are  entered  and  the  filter  goes  through  sufficient  iteration  to  converge. 
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IV.  FUTURE  DEVELOPMENT 


A.  SUGGESTIONS  FOR  PROGRAM  IMPROVEMENT 

LosCon’s  interface  is  a  simple  Microsoft  Excel  spreadsheet  using  the  graphical 
plug-in  to  generate  visual  representations  of  the  mathematical  program  results.  These 
built-in  graphical  capabilities  do  not  use  standard  Navy  icons.  An  updated  version  could 
use  a  different  interface  to  allow  for  surface  ships,  submarines,  and  aircraft  to  be 
displayed  in  standard  icons.  While  creating  an  independent  interface  program  may  be 
simpler,  keeping  the  program  essentially  an  Excel  file  allows  it  to  be  easily  burned  on 
CDs  and  used  as  a  source  of  additional  information  during  tactical  engagements  and 
exercises. 

In  some  ways,  LosCon  may  be  too  automated.  Making  it  possible  for  the  operator 
to  change  some  built-in  variables  from  their  default  settings  could  make  LosCon  more 
useful.  Changing  the  root  mean  square  velocity  to  correspond  to  specific  submarine 
types  would  increase  the  usefulness  of  LosCon.  Tables  with  cruising  speeds  for  each 
type  of  submarine  could  be  added  to  the  program,  so  an  operator  could  input  the 
submarine  type  and  LosCon  would  call  the  appropriate  root  mean  square  velocity  for  use 
in  calculations.  The  default  long-tenn  average  velocity  for  the  enemy  submarine  could 
also  be  made  variable  to  allow  LosCon  to  be  used  for  tracking  transiting  submarines. 

LosCon  could  also  be  made  easier  to  use  if  parts  of  it  were  more  automated.  The 
iteration  process  to  find  an  initial  estimate  that  converges  could  be  built  into  the  code  so 
that  it  was  accomplished  with  one  click  rather  than  many.  The  re-initialization  of  the 
initial  estimate  with  every  new  contact  report  could  also  be  built  in.  Ideally,  one  button 
press  would  start  the  iteration  process  to  generate  a  converging  initial  estimate  based  on 
the  first  few  contact  reports,  and  then  apply  all  the  additional  measurements  to  that 
estimate  to  generate  an  AOU  at  the  time  of  the  last  contact  report. 

For  LosCon  to  be  used  optimally,  the  single  standard  deviation  error  needs  to  be 
included  in  the  program  input.  Standard  contact  report  formats  used  on  surface  ships  do 
not  require  inclusion  of  the  standard  deviations  of  bearing  error  or  range  error.  The  errors 
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are  not  so  random  that  standard  deviations  are  unattainable.  Towed  arrays,  for  example, 
have  varying  errors  based  on  the  directionality.  The  exact  numbers  would  vary  to  some 
small  extent  but  the  test  values  generated  when  the  system  was  first  installed  could  be 
used.  Error,  even  when  it  is  reported,  is  not  in  a  standard  format  such  as  a  standard 
deviation  or  variance.  It  is  important  that  the  single  standard  deviation  values  inputted  in 
LosCon  be  more  than  operator  estimates.  If  standardization  of  error  reports  were 
sufficiently  emphasized  for  good  records  to  be  kept  of  current  and  developing  systems, 
single  standard  deviation  bearing  errors  in  degrees  based  on  source  and  directionality 
could  be  tabulated  and  included  in  the  next  version  of  LosCon. 

B.  TACTICAL  DECISION  AID  DEVELOPMENT  OPTIONS 

The  usefulness  of  LosCon  is  limited,  because  it  does  not  deal  with  ocean 
acoustics.  A  number  of  questions  involving  both  oceanographic  and  probabilistic  data 
remain  to  be  answered.  Where  inside  the  LosCon-generated  AOU  are  operators  able  to 
pick  up  sonar  signals,  for  instance?  Where  should  surface  platforms  reposition 
themselves  to  form  sonar  nets  with  onboard  sensors  and  towed?  Where  should  maritime 
patrol  craft  drop  sonobuoy  fields  to  optimize  the  chances  of  reacquiring  the  lost  contact? 
An  ideal  ASW  tactical  decision  aid  package  would  include  both  oceanographic  and 
probabilistic  data. 

The  properties  of  the  ocean  which  affect  sound  propagation  and  thus  sonar 
effectiveness  have  been  well  studied  and  well  modeled  by  a  number  of  software 
packages — perhaps  the  most  well  known  is  IMAT.24  These  programs  deal  in  the 
oceanographic  issues  for  optimizing  the  volume  of  water  with  can  be  observed  through 
sonar.  These  are  traditionally  extremely  data-intensive  computations  that  require 
dedicated  hardware.  It  is  still  a  data-intensive  process,  but  mathematical  and  coding 
efforts  have  slimmed  the  program  enough  to  create  a  useable  shipboard  oceanographic 
modeling  program  in  PC  IMAT. 


24  Czech,  Carl.  U.S.  Navy  Publication,  NPRDC-TR-98-2  The  Interactive  Multisensor  Analysis 
Training  (IMAT)  System:  an  Evaluation  of  Airborne  Acoustic  Mission  Course,  (San  Diego,  CA:  GPO, 
1998). 
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The  nature  of  Kalman  filtering  in  which  data  is  repeatedly  overwritten  makes 
memory  demands  of  generating  the  AOU  for  a  contact  extremely  small.  LosCon  fits  on  a 
floppy  disk  even  with  contact  reports.  If  LosCon  were  compiled  with  a  new  version  of 
PC  IMAT  or  a  similar  oceanographic  program,  LosCon  would  add  negligible  amounts  to 
system  memory  requirements  and  file  size.  The  value  of  the  final  product,  however, 
would  be  greatly  increased. 

The  ASW  commander  could  look  at  a  single  screen  displaying  the  convergence 
zone  rings  from  combatant  platfonns  and  offboard  sensors  along  with  the  enemy 
submarine’s  area  of  uncertainty.  The  commander  could  then  overlay  sonobuoy  patterns 
and  reposition  assets  over  the  AOU  to  maximize  the  likelihood  of  reacquiring  contact. 

Such  an  ASW  package  could  easily  include  the  following  displays: 

•  Display  the  effects  of  oceanographic  properties  on  sonar  equipment.  (IMAT)  25 

•  Display  areas  of  uncertainty  for  multiple  contact  reports.  (LosCon) 

•  Display  AOUs  at  a  future  time  t  if  no  new  contact  reports  are  made.  (LosCon) 

Once  probabilistic  target  location  and  oceanographic  properties  are  combined,  a 
myriad  of  extremely  useful  options  become  possible.  Given  a  time  for  the  drop  and  the 
number  of  sonobuoys  available,  one  could  generate  the  optimal  sensor  placement.  That 
sonobuoys  could  only  be  dropped  in  a  fixed  number  of  lines  could  be  included  in  the 
constraints  along  with  an  expected  variance  to  the  exact  placement  plan.  Optimal  sensor 
placement  could  be  generated,  given  additional  limits  such  as  not  allowing  a  carrier  to 
change  course  or  speed  during  a  fixed  period  due  to  flight  operations.  Loss  of  contact 
itself  could  become  a  tracking  tool  if  the  combination  of  a  submarine’s  probabilistic 
location  and  the  ocean  environment  allowed  the  war-fighter  to  predict  that  the  AOU 
would  pass  out  of  the  convergence  zone. 

LosCon  is  a  useful  and  proven  program  which  can  help  war-fighters  regain 
tactical  control  in  a  loss  of  contact  situation.  Persistent  detection  and  cueing,  one  of  the 
fundamental  principles  of  Chief  of  Naval  Operations  Admiral  Clark’s  ASW  Concept  of 
Operations,  depend  on  utilizing  contact  reports  from  a  network  of  combatant  platforms 

25  Czech,  Carl.  U.S.  Navy  Publication,  NPRDC-TR-98-2  The  Interactive  Multisensor  Analysis 
Training  (IMAT)  System:  an  Evaluation  of  Airborne  Acoustic  Mission  Course,  (San  Diego,  CA:  GPO, 
1998). 
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and  offboard  sensors.26  At  present,  LosCon  can  integrate  past  ASW  contact  information 
to  provide  war-fighters  with  real-time  probabilistic  location  of  target  submarines.  In  the 
future,  LosCon  could  be  part  of  a  tactical  decision  aid  that  makes  sea  shield  a  reality. 


26  Clark,  Vern.  Anti-Submarine  Warfare  Concept  of  Operations  (Draft  26  April  2004). 
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APPENDIX  A:  USER’S  MANUAL 


A.  SELECTING  A  MAP  ORIGIN 

A  latitude  and  longitude  near  the  area  of  operations  is  ideal  for  a  base  position. 
This  selection  can  be  changed  as  desired  by  the  operator,  but  all  the  target  positions  must 
be  updated  after  doing  so.  In  the  graph  of  target  positions,  the  reference  latitude  and 
longitude  serve  as  the  map  origin,  so  select  a  set  of  coordinates  in  the  southwest  corner  of 
the  area  of  operations  to  keep  all  information  in  the  first  quadrant.  Figure  7  shows  where 
to  enter  the  reference  coordinates. 
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Target3 
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80.0000 

lat-y 
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6.3773 

80.0000 
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Update] 

Update2 

|  Update3 

Latitude  (deg) 

35.0000 

35.0000 

36.0000 

Latitude  (min) 

5.5171 

6.3773 

20.0000 

Longitude  (deg) 

119.0000 

113.0000 

118.0000 

Longitude  (min) 

50.0758 

42.6998 

22.3380 

refjat  (deg) 

35 

ref  Ion  (deg  W) 

120 

Figure  7:  The  map  origin  is  entered  as  reference  latitude  and  reference 

longitude  in  the  master  spreadsheet. 


B.  ENTERING  CONTACT  REPORTS 

Contact  reports  are  entered  in  the  master  spreadsheet  in  two  formats.  A  position 
can  be  entered  directly  as  a  latitude  and  longitude  of  the  target  with  associated  circularly- 
symmetric  standard  deviation  error  in  nautical  miles,  or  a  ship  can  enter  its  own  latitude 
and  longitude  combined  with  the  bearing  to  the  contact  and  standard  deviation  bearing 
error  as  can  be  seen  in  Figure  8. 
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Source 
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35 

45 
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2 

45 
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1 

0 

i 
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3 

35 

0 
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0 

0 

0 
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1 

2 

2 

0.0024 

Figure  8:  Contact  reports  are  entered  into  the  master  spreadsheet. 


To  enter  a  bearing  contact  report,  enter  the  latitude  and  longitude  of  the  ship  that 
took  the  bearing  in  columns  three  through  six.  Do  not  enter  anything  in  columns  seven 
and  eight.  The  north  south  and  east-west  coordinates  will  automatically  be  calculated  in 
columns  seven  and  eight,  respectively.  Enter  the  bearing  to  the  contact  in  column  nine  in 
degrees.  Enter  the  single  standard  deviation  bearing  error  in  column  ten.  In  column 
eleven,  enter  on  which  target  the  contact  report  is  believed  to  be.  Enter  the  time  of  the 
contact  report  in  column  twelve  in  terms  of  hours  since  the  beginning  of  the  tracking 
exercise.  Enter  a  “1”  in  column  thirteen  indicating  to  LosCon  that  the  contact  report 
entered  on  that  row  should  be  included  as  a  bearing  rather  than  a  position  in  the  next 
update  of  the  target’s  location.  LosCon  calculates  the  dimensionless  shock  in  column 
fourteen. 

To  enter  a  position  contact  report,  enter  the  latitude  and  longitude  of  the  position 
in  columns  three  through  six.  As  before,  the  north  south  and  east-west  coordinates  will 
automatically  be  calculated  in  columns  seven  and  eight,  respectively.  Leave  the  bearing, 
column  nine,  blank.  Enter  the  single  standard  deviation  circularly  symmetric  error  for  the 
position  in  nautical  miles  in  column  ten.  In  column  eleven,  enter  on  which  target  the 
contact  report  is  believed  to  be.  Enter  the  time  of  the  contact  report  in  column  twelve, 
and  enter  a  “2”  in  column  thirteen  indicating  to  LosCon  that  the  contact  report  entered  on 
that  row  should  be  included  as  a  position  in  the  next  update  of  the  target’s  location. 

Again,  LosCon  calculates  the  dimensionless  shock  in  column  fourteen. 

Dimensionless  shock  should  be  used  as  a  warning  signal.  The  dimensionless 
shock  can  be  thought  of  as  how  much  the  program  is  surprised  by  the  new  information 
the  report  provides  on  the  submarine’s  location.  Large  dimensionless  shock  does  not 
mean  that  a  contact  report  is  necessarily  false,  but  a  false  contact  report  will  probably 
have  a  large  dimensionless  shock.  If  the  past  few  contact  reports  were  false  or  had  large 
errors  in  the  same  direction,  the  next  true  contact  report  would  have  a  large  dimensionless 
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shock.  These  are  divergence  issues.  See  “E:  Divergence  and  Convergence”  for  more 
information. 


C.  ESTIMATING  INITIAL  TARGET  LOCATION 

The  operator  must  estimate  an  initial  position  for  each  submarine.  The  position 
must  be  within  a  few  standard  deviations  of  the  first  few  contact  reports  for  the 
submarine.  An  estimate  too  far  from  the  earliest  contact  reports  may  cause  the  program 
to  return  an  infinitely  large  AOU  (Area  Of  Uncertainty)  as  the  program  is  iterated.  This 
is  referred  to  as  divergence.  The  operator  can  change  the  initial  position  for  each 
submarine  at  any  time,  but  after  selecting  an  initial  position  that  is  feasible  compared  to 
the  first  few  contact  reports,  there  is  no  reason  to  use  a  different  initial  position. 

The  estimated  location  and  velocity  of  the  target  submarine  is  in  the  four  by  one 
matrix,  highlighted  in  Figure  9.  From  top  to  bottom,  the  cells  are  east-west  position  in 
nautical  miles  relative  to  the  selected  origin,  north-south  position  in  nautical  miles,  east- 
west  velocity  component  in  knots,  and  north-south  velocity  component  in  knots.  The 
four  by  four  matrix,  also  highlighted,  is  the  covariance  matrix  ordered  in  the  same  way  as 
the  location  and  velocity  matrix.  The  diagonals  of  the  covariance  matrix  are  variances. 
The  latitude  and  longitude  coordinates  for  the  target  location  are  automatically  calculated 
from  the  values  entered  for  the  east- west  position  and  north-south  position.  If  the 
reference  latitude  and  longitude  were  changed  on  the  master  spreadsheet,  the  latitude  and 
longitude  that  is  five  nautical  miles  east  of  the  origin  and  five  nautical  miles  north  of  the 
origin  would  change  to  correspond.  The  “Update  to”  button  is  used  for  generating  an 
AOU  at  some  time  after  the  last  contact  report.  See  “D:  Estimating  an  AOU.” 
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Figure  9:  An  initial  estimate  for  the  location,  speed,  and  accuracy  of  that 

estimate  is  entered  in  spreadsheet,  “Targetl.” 

An  initial  estimate  for  the  location  of  each  submarine  must  be  entered  in  LosCon. 
This  value,  if  reasonable,  will  be  changed  to  the  optimal  initial  estimate  in  the  course  of 
iteration.  The  next  ten  figures  will  show  how  to  iterate  the  extended  Kalman  filter  until  it 
converges  on  an  optimal  initial  estimate.  As  demonstrated  in  Figure  10,  an  initial 
estimate  for  the  location  of  each  target  submarine  should  be  entered  in  the  two  cells 
immediately  below  “Working  Estimate.”  These  are  east-west  and  north-south 
coordinates  in  nautical  miles.  The  origin  is  the  one  selected  on  the  “Master”  spreadsheet. 
See  “A:  Selecting  a  Map  Origin.”  The  latitude  and  longitude,  which  correspond  to  the 
estimated  location,  will  appear  in  degrees  and  minutes  next  to  the  two  entries. 
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Figure  10:  An  initial  estimate  is  entered  in  spreadsheet  “Targetl”  as  (x  = 
0  nm,  y  =  0  nm)  and  a  velocity  of  zero.  The  position  corresponds  to  (35  °N,  120  °W), 
the  selected  map  origin. 


Use  the  first  few  bearing  or  position  reports  to  find  the  optimal  initial  estimate  for 
the  location  of  a  submarine.  Figure  1 1  shows  two  sample  bearing  contact  reports.  Based 
on  these  two  measurements,  the  initial  position  estimate  will  be  iterated  until  it  converges 
or  diverges  after  pressing  the  “Update  1”  button  repeatedly.  Figures  12  through  15  show 
the  values  of  the  position  estimates  based  on  the  contact  reports  from  Figure  11.  That 
position  estimate  is  iterated  until  it  converges.  See  “D:  Divergence  and  Convergence.” 
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Figure  11:  The  two  bearing  contact  reports  are  entered  for  “Target  #1” 
on  the  “Master”  spreadsheet. 
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Figure  12:  The  estimated  location  for  “Target  #1”  on  the  “Master” 
spreadsheet  after  entering  an  estimate  is  the  same  as  was  entered  on  the  “Targetl” 
spreadsheet. 
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Figure  13:  The  estimated  location  for  “Target  #1”  changes  after  the 
“Updatel”  button  is  pressed  once. 
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Figure  14:  The  estimated  location  for  “Target  #1”  continues  to  change  the 
second  time  the  “Updatel”  button  is  pressed. 
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Figure  15:  After  the  “Updatel”  button  is  pressed  nine  times  the  estimated 
target  location  converges  at  (35°  5.6’  N,  119°  50’  W). 


After  an  optimal  initial  estimate  for  the  location  of  a  submarine  is  found,  that 
value  needs  to  be  saved  for  future  use.  Every  time  a  new  contact  report  on  that 
submarine  is  entered,  the  optimal  initial  estimate  should  be  re-entered  as  well.  Figure  16 
shows  the  two  position  entries  from  the  optimal  initial  estimate  being  saved. 
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Figure  16:  The  optimized  initial  estimate  in  spreadsheet  “Targetl”  is 
copied  from  the  “Working  Estimate”  cells  to  the  “Initial  Location  Estimate”  cells 
for  later  use. 


A  new  contact  report  for  “Target  #1”  is  shown  in  Figure  17.  The  estimate  for  the 
location  of  the  submarine  no  longer  needs  to  be  iterated,  so  the  “Update  1”  button  should 
be  pressed  only  once.  Figures  18  and  19  show  how  that  new  contact  report  changes  the 
estimate  location  at  the  time  of  the  third  contact  on  the  “Master”  and  “Targetl” 


spreadsheets. 
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Figure  17:  A  new  “Target  #1”  contact  report  for  a  time  two  hours  later  is 
entered. 
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Figure  18:  The  estimated  location  of  the  submarine  is  shown  on 
spreadsheet  “Master”  after  entering  a  third  contact  report  and  pressing  “Updatel” 
once. 
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Figure  19:  The  estimated  location  of  the  submarine  is  also  shown  on 
spreadsheet  “Targetl”  after  entering  a  third  contact  report  and  pressing  “Updatel” 
once. 


The  time  in  hours  shown  in  the  two  cells  above  the  “Update  to”  button  is  the  time 
of  the  last  contact  report  and  is  the  time  for  which  the  “Working  Estimate”  of  submarine 
location  is  valid.  The  “Working  Estimate”  may  not  be  a  good  estimate  for  time  zero,  so 
as  new  contacts  are  entered  on  the  “Master”’  spreadsheet,  the  values  saved  in  “Initial 
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Location  Estimate”  should  be  entered  in  the  two  cells  in  the  box  immediately  below 
“Working  Estimate”  before  pressing  the  “Update”  button. 


D.  ESTIMATING  AN  AOU 

On  each  target-specific  spreadsheet,  the  highlighted  cell  above  the  button, 
“Update  to”  contains  the  time  of  the  AOU  currently  being  displayed.  The  AOU  at  the 
time  of  the  last  contact  report  is  automatically  calculated  as  shown  in  Figure  20.  The 
center  of  the  AOU  is  the  position  shown  under  “Working  Estimate.”  The  orientation  in 
degrees  is  the  bearing  of  the  AOU’s  major  axis.  The  value  under  “major(nm)”  is  the 
major  axis’s  length  in  nautical  miles.  Likewise,  the  value  under  “minor(nm)”  is  the 
minor  axis’s  length  in  nautical  miles.  The  minor  axis  is  perpendicular  to  the  major  axis, 
and  thus  it  has  an  orientation  ninety  degrees  off  the  orientation  listed. 

To  generate  an  AOU  for  a  time  after  the  time  of  the  last  contact  report,  do  not 
directly  change  the  number  in  the  cell  which  displays  the  time  of  the  last  contact  report. 
That  will  generate  an  inaccurate  result,  because  the  other  numbers  used  in  the  calculation 
will  not  have  been  updated.  Entering  the  time  desired  in  hours  in  the  cell  above  the 
“Update  to”  button,  selecting  a  different  cell,  and  then  clicking  on  “Update  to”  will  move 
the  number  entered  up  one  row  and  recalculate  that  submarine’s  AOU  for  the  desired 
time.  Figure  20  and  Figure  21  show  an  AOU  being  updated  to  a  future  time. 
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Figure  20:  Changing  the  time  generates  an  updated  AOU. 
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Figure  21:  The  matrix  is  recalculated  after  the  “Update  to”  button  is 
clicked. 


E.  DIVERGENCE  AND  CONVERGENCE 

When  LosCon  is  being  used  optimally,  the  iteration  of  the  filter  (repeatedly 
pressing  the  “Update”  button)  will  cause  convergence.  This  means  that  with  each 
iteration,  the  amount  of  change  to  the  center  point  of  the  AOU  will  decrease.  Iteration  is 
desirable,  because  it  increases  the  accuracy  of  the  program  by  minimizing  the  effect  of 
the  initial  estimate  of  the  target’s  position.27  The  goal  is  to  generate  an  “Initial  Location 
Estimate”  by  iterating  the  filter  with  the  first  few  contact  reports.  Those  values  should  be 
used  for  the  initial  estimate  every  time  a  new  contact  report  is  entered  in  LosCon.  See 
“C:  Estimating  Initial  Target  Locations”  for  details  on  when  and  how  to  iterate  the 
estimated  target  location. 

Divergence  refers  to  the  event  where  the  extended  Kalman  filter  fails  to  correct 
errors  in  the  appropriate  direction  and  instead  compounds  them  by  moving  the  target 
location  estimate  farther  and  farther  from  the  true  location  as  the  filter  is  iterated. 
Divergence  can  be  caused  by  a  poor  initial  target  estimate  or  a  series  of  false  contact 
reports  that  are  inputted  as  true.  When  iteration  does  not  result  in  decreasing  changes  in 
the  location  of  the  center  of  the  AOU,  divergence  has  occurred. 

If  divergence  occurs,  a  new  initial  estimate  for  the  location  and  velocity  of  the 
target  must  be  entered.  LosCon’s  working  estimate  for  the  location  of  the  submarine  is  a 

27  Kerr,  Thomas  H..  “Streamlining  Measurement  Iteration  for  EKF  Target  Tracking,”  IEEE 
Transactions  on  Aerospace  and  Electronic  Systems  27,  no.  2  (March  1991):  408-421. 
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product  of  the  last  update.  If  the  program  has  just  diverged,  that  working  estimate  is 
probably  an  extremely  bad  initial  estimate.  Changing  the  initial  estimate  alone  may  be 
sufficient  to  reacquire  a  convergence. 

If  entering  a  new  initial  estimate  alone  is  not  sufficient,  again  enter  a  feasible 
initial  estimate  and  select  a  different  set  of  early  contact  reports  to  use  in  the  next  attempt 
to  find  a  converging  estimate  for  the  location  of  the  target  submarine  at  time  zero.  Recall 
that  LosCon  allows  reports  to  be  included  and  left  out  with  a  simple  toggle  switch  at  any 
time.  Trying  a  variety  of  sets  of  reports  and  reentering  a  feasible  initial  estimate  each 
time,  should  result  in  converging  output. 

F:  OPERATION  SUMMARY 

LosCon’ s  iterated  extended  Kalman  filter,  which  generates  the  center  coordinate 
of  the  AOU,  works  in  Cartesian  coordinates.  The  operator  must  enter  a  reference  latitude 
and  longitude  point  for  the  program  to  minimize  the  error  associated  with  converting 
approximately  spherical  coordinates  to  locally-flat  space.  LosCon  automatically  converts 
the  latitude  and  longitude  coordinates  entered  by  the  operator  in  contact  reports  to  east- 
west  and  north-south  coordinates  in  nautical  miles  referenced  to  the  origin  chosen  by  the 
operator.  The  outputs  are  presented  in  Cartesian  coordinates  and  in  latitude  and 
longitude.  See  “A:  Selecting  a  Map  Origin.” 

Contact  reports  can  be  bearings  or  positions.  The  contact  report  database  remains 
available  to  the  operator  to  change  and  or  update  as  desired.  Contact  reports  ideally 
consist  of  the  raw  data:  the  latitude  and  longitude  of  the  observer  combined  with  a 
bearing  on  the  target  and  an  estimated  standard  deviation  bearing  error  at  a  given  time.  If 
the  contact  report  is  instead  a  position  derived  from  TMA  (Target  Motion  Analysis), 
LosCon  will  accept  a  target  latitude  and  longitude  combined  with  that  position’s  standard 
deviation  circularly  symmetric  error  in  nautical  miles  and  the  time  of  the  report.  See  “B: 
Entering  Contact  Reports.” 

On  the  subsidiary  spreadsheets  assigned  to  each  potential  target,  the  operator  must  enter 
an  initial  estimate  for  the  enemy  submarines’  locations  at  time  zero.  A  good  estimate  is  a 
latitude  and  longitude  near  the  search  assets  but  not  exactly  the  same.  If  a  submarine  is 

found  it  will  be  relatively  near  the  searchers,  so  that  initial  estimate  is  logical  and  unlikely 
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to  make  the  search  problem  diverge.  If  the  first  estimate  results  in  divergence,  the 
operator  can  simply  try  a  new  one  a  few  seconds  later.  The  initial  estimate  includes  two 
velocity  components  as  well  as  location;  these  velocity  estimates  do  not  affect 
convergence  and  should  usually  be  set  to  zero.  See  “C:  Estimating  Initial  Target 
Location”  and  “E:  Divergence  and  Convergence.” 

The  AOUs  generated  on  the  target  specific  spreadsheets  correspond  to  each 
individual  target.  The  centers  of  the  AOUs  are  updated  on  the  master  spreadsheet  as 
well.  AOUs  can  also  be  calculated  for  future  times.  See  “D:  Estimating  an  AOU.” 
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APPENDIX  B:  VISUAL  BASIC  CODE 


A.  MASTER  CALCULATIONS 

1.  Description  of  Code 

LosCon’s  code  is  after  the  MTST  (Maneuvering  Target  Statistical 
Tracker)  program  by  Professor  Alan  Washburn  of  the  Naval  Postgraduate  School. 
Daniel  H. Wagner  designed  the  first  MTST  program  in  1980.  A  copy  of 
LosCon.xls  can  be  obtained  from  Prof.  Washburn  (email:  awashburn@nps.edu). 

The  following  code  has  been  modified  to  include  three  distinct  targets, 
false  contact  reports,  position  and  bearing  target  types,  and  latitude  and  longitude 
inputs  vice  Cartesian  coordinates.  Changing  the  default  single  standard  deviation 
for  the  root  mean  square  of  the  target  submarines’  velocities  and  changing  the 
default  long-term  average  velocity  of  the  submarines  can  be  done  in  Kmove 
subroutine. 


2.  Initialization  and  Cartesian  Conversion 


a.  Initialization  and  Global  Variables 

Option  Base  1 
Option  Explicit 
Dim  Sigma,  Mu 

b.  Latitude  to  Cartesian  Coordinate  Conversion 

Function  Latdeg(y,  latref,  lonref) 

Dim  latitude 
latitude  =  latref  +  y  /  60 
Latdeg  =  Int(latitude) 

End  Function 

Function  Latmin(y,  latref,  lonref) 

Dim  latitude 

latitude  =  latref  +  y  /  60 

Fatmin  =  (latitude  -  Int(latitude))  *  60 
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End  Function 


Public  Function  yFat(Fatdeg,  Fatmin,  latref,  lonref) 

Dim  latitude 

latitude  =  Fatdeg  +  Fatmin  /  60yFat  =  (latitude  -  latref)  *  60 
End  Function 

c.  Longitude  to  Cartesian  Coordinate  Conversion 

Function  Fondeg(x,  latref,  lonref) 

Dim  longitude,  fac 
fac  =  60  *  Cos(Atn(l)  *  latref  /  45) 
longitude  =  lonref  -  x  /  fac 
Fondeg  =  Int(longitude) 

End  Function 

Function  Fonmin(x,  latref,  lonref) 

Dim  longitude,  fac 

fac  =  60  *  Cos(Atn(l)  *  latref  /  45) 

longitude  =  lonref  -  x  /  fac 

Fonmin  =  (longitude  -  Int(longitude))  *  60 

End  Function 

Public  Function  xFon(Fondeg,  Fonmin,  latref,  lonref) 

'longitudes  are  in  degrees  west 

’xFon  is  positive  if  longitude  is  east  of  the  reference. 

Dim  longitude,  fac 
fac  =  60  *  Cos(Atn(l)  *  latref  /  45) 
longitude  =  Fondeg  +  Fonmin  /  60 
xFon  =  (lonref  -  longitude)  *  fac 
End  Function 


3.  Updating  Target  Submarines 

Public  Sub  Update(target) 

Dim  Minputs,  Mu,  Sigma 
Dim  nrows&,  ncols&,  sht$ 

Select  Case  target 
Case  1 

sht  =  "  Target  1" 

With  Worksheets(sht) 

Sigma  =  .Range("MATl").  Value 
Mu  =  .Range("muin").  Value 
End  With 
Case  2 

sht  =  "Target2" 

With  Worksheets(sht) 
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Sigma  =  .Range("MATl").  Value 
Mu  =  .Range("muin").  Value 
End  With 
Case  3 

sht  =  "Target3" 

With  Worksheets(sht) 

Sigma  =  .Range("MATl").  Value 
Mu  =  .Range("muin").  Value 
End  With 
Case  Else 
Exit  Sub 
End  Select 

With  Range("MasterIn") 

Minputs  =  .Value 
mows  =  .Rows. Count 
ncols  =  .Columns.  Count 
With  .Interior 
.Colorlndex  =  35 
.Pattern  =  xlSolid 
End  With 
End  With 
Dim  i&,  j&,  k&,  t# 

For  i  =  1  To  4 
For  j  =  1  To  4 
Sigma(i,  j)  =  0 
Next  j 
Next  i 

Sigma(  1 ,  1)  =  100000:  Sigma(2,  2)  =  100000 
Sigma(3,  3)  =  100:  Sigma(4,  4)  =  100 
t  =  0 

For  k  =  1  To  nrows 
'If  include  doesn’t  =  0 
If  Minputs(k,  12)  >  0  Then 

'If  this  is  the  target  you  are  processing 
If  Minputs(k,  10)  =  target  Then 
Call  KMove(Mu,  Sigma,  Minputs(k,  11)  -  t) 
t  =  Minputs(k,  1 1) 

If  Minputs(k,  12)  =  1  Then 

Call  KMeas(Mu,  Sigma,  Minputs(k,  7),  Minputs(k,  6),  Minputs(k,  8), 
Minputs(k,  9),  Minputs(k,  13)) 

Else 

Call  KMeasPos(Mu,  Sigma,  Minputs(k,  7),  Minputs(k,  6),  Minputs(k,  9), 
Minputs(k,  13)) 

End  If 

Range("MasterIn")(k,  13).Value  =  Minputs(k,  13) 

End  If 
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'For  not  included  rows 
Else 

Range("MasterIn")(k,  13).Value  =  "" 

End  If 
Next  k 

With  Worksheets(sht) 

.Range("MATl").  Value  =  Sigma 
.Range("muin").  Value  =  Mu 
.Range("CurTime")(l,  1). Value  =  t 
End  With 

’Worksheets("Master").Range(sht).Value  =  Mu 
End  Sub 

4.  Bearing  Measurements 

Sub  KMeas(Mu,  Sigma,  x,  y,  theta,  sig,  dshock) 

Dim  dx#,  dy#,  ran#,  ang#,  fac#,  H(  1 ,  4)  As  Double 
Dim  Shock#,  r#,  el,  e2,  denom,  KGain(4,  1)  As  Double 
Dim  lat#,  Ion#,  latmu#,  lonmu#,  latref#,  lonref# 

Dim  beta#,  conv# 
conv  =  Atn(l)  /  45 

latref  =  CDbl(Worksheets("master").Range("p  13").  Value) 
lonref  =  Worksheets("master").Range("pl4").Value 
conv  =  1  /  (60  *  Cos(latref  *  conv)) 
lat  =  (latref  +  y  /  60) 

Ion  =  (lonref  -  conv  *  x) 
latmu  =  (latref  +  Mu(2,  1)  /  60) 
lonmu  =  (lonref  -  conv  *  Mu(l,  1)) 
dx  =  Mu(l,  1)  -  x:  dy  =  Mu(2,  1)  -  y 
ran  =  dx  *  dx  +  dy  *  dy 
If  ran  =  0  Then  Exit  Sub 
r  =  (sig  *  conv)  A  2 

'  following  two  statements  replaced  by  bearing()  call  4/21/04 
'ang  =  Atn2(dy,  dx) 

'Shock  =  theta  *  conv  -  ang 

ang  =  bearing(lat,  lonmu,  latmu,  Ion)  'note  Ion  reversal 
Shock  =  (theta  -  ang)  *  conv 
If  Shock  >  4  *  Atn(l)  Then 
Shock  =  Shock  -  8  *  conv 
Elself  Shock  <  -4  *  Atn(l)  Then 
Shock  =  Shock  +  8  *  conv 
End  If 

H(l,  1)  =  dy  /  ran:  H(l,  2)  =  -dx  /  ran 
H(l,  3)  =  0:  H(l,  4)  =  0 

el  =  Application.  WorksheetFunction.Transpose(H) 
e2  =  Application. WorksheetFunction.MMult(Sigma,  el) 
denom  =  Application.WorksheetFunction.MMult(H,  e2) 
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dshock  =  Shock  *  Shock  /  (denom(l)  +  r) 

If  dshock  >  1000  Then 

Exit  Sub  ’but  dshock  is  still  returned 
End  If 

'update  mu  and  sigma 
fac  =  (denom(l)  +  r)  A  (-1) 

'  can't  figure  out  how  to  make  a  function  return  an  array,  so  use  subs 
’  last  argument  in  subs  is  the  output 
Call  matscal(e2,  fac,  KGain) 

Call  matscal(KGain,  Shock,  el) 

Call  matadd(Mu,  el,  Mu) 

el  =  Application. WorksheetFunction.MMult(KGain,  H) 
e2  =  Application.WorksheetFunction.MMult(el,  Sigma) 

Call  matdif(Sigma,  e2,  Sigma) 

'Debug.Print  dx;  dy;  ran 
End  Sub 

5,  Position  Measurements 

Sub  KMeasPos(Mu,  Sigma,  x,  y,  sig,  dshock) 

'input  is  a  target  location,  not  a  bearing 

'x  and  y  are  the  measured  target  location 

'sig  is  the  standard  deviation  of  x  and  of  y,  in  compatible  units 

Dim  q#,  H(2,  4)  As  Double 

Dim  Shock(2,  1)  As  Double,  el,  e2,  denom,  KGain 

Shock(l,  1)  =  x  -  Mu(l,  1):  Shock(2,  1)  =  y  -  Mu(2,  1) 

q  =  sig  *  sig 

H(l,  1)  =  1:  H(l,  2)  =  0:  H(l,  3)  =  0:  H(l,  4)  =  0 
H(2,  1)  =  0:  H(2,  2)  =  1:  H(2,  3)  =  0:  H(2,  4)  =  0 
el  =  Application.  WorksheetFunction.Transpose(H) 
e2  =  Application. WorksheetFunction.MMult(Sigma,  el) 
denom  =  Application.WorksheetFunction.MMult(H,  e2) 
denom(l,  1)  =  denom(l,  1)  +  q:  denom(2,  2)  =  denom(2,  2)  +  q 
denom  =  Application.WorksheetFunction.MInverse(denom) 
el  =  Application. WorksheetFunction.MMult(denom,  Shock) 
dshock  =  Shock(l,  1)  *  el(l,  1)  +  Shock(2,  1)  *  el(2,  1) 

If  dshock  >  1000  Then 

Exit  Sub  'but  dshock  is  still  returned 
End  If 

KGain  =  Application.WorksheetFunction.MMult(e2,  denom) 
el  =  Application. WorksheetFunction.MMult(KGain,  Shock) 

'  can't  figure  out  how  to  make  a  function  return  an  array,  so  use  subs 

'  last  argument  in  subs  is  the  output 

Call  matadd(Mu,  el,  Mu)  'Mu  is  now  updated 

el  =  Application. WorksheetFunction.MMult(KGain,  H) 

e2  =  Application.WorksheetFunction.MMult(el,  Sigma) 

Call  matdif(Sigma,  e2,  Sigma)  'Sigma  is  now  updated 
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End  Sub 


6.  MTST  Movement  Subroutine 

Sub  KMove(Mu,  Sigma,  t  As  Double) 

Dim  s#,  beta#,  bl#,  b2#,  c#,  cl#,  c2#,  c3# 
s  =  10:  beta  =  2 

's  is  rms  speed  (kt),  1/beta  is  the  relaxation  time  (hr) 

’To  change  s  the  root  mean  square  speed  of  the  target  submarines 
'to  potentially  different  rms  speeds  for  each,  replace  the  line  s=10  with 
'Dim  si#,  s2#,  s3# 

'si  =  Worksheets("Targetl").Range("name  of  cell"). Value 
's2  =  Worksheets("Target2").Range("name  of  cell").  Value 
's3  =  Worksheets("Target3").Range("name  of  cell").  Value 
'where  name  of  cell  is  a2  or  h32  for  example. 

'Be  sure  to  keep  the  quotation  marks  when  choosing  a  cell  in  which  to 
'input  the  rms  speeds. 

Dim  matin,  matout,  el,  e2,  e3,  e4 
Dim  i&,  j&,  k& 

'Dim  el (4,  4)  As  Double,  e2(4,  4)  As  Double,  e3(4,  4)  As  Double,  e4(4,  4)  As 
Double 

Dim  Phi(4,  4)  As  Double,  q(4,  4)  As  Double 
'matin  =  mat.  Value 

'matout  =  Application.WorksheetFunction.MInverse(matin) 

For  i  =  1  To  4 
For  j  =  1  To  4 
Phi(i,j)  =  0:  q(i,j)  =  0 
Next  j 
Next  i 

If  beta  >  0  Then 

If  t  *  beta  <  700  Then 
b2  =  Exp(-t  *  beta) 

Else 
b2  =  0 
End  If 

b  1  =  ( 1  -  b2)  /  beta 
c  =  0.5  *  s  *  s 

cl  =  c*(t  +  t-(3-4*b2  +  b2*  b2)  /  beta)  /  beta 
c2  =  c  *  beta  *  b  1  *  b  1 
c3  =  c  *  ( 1  -  b2  *  b2) 

Else 
b2  =  1 
bl  =t 
cl  =  0 
c2  =  0 
c3  =  0 
End  If 
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Phi(  1 ,  1)  =  1:  Phi(2,  2)  =  1 
Phi(3,  3)  =  b2:  Phi(4,  4)  =  b2 
Phi(l,  3)  =  bl:  Phi(2,  4)  =  bl 
q(l,  1)  =  cl:  q(2,  2)  =  cl 
q(3,  3)  =  c3:  q(4,  4)  =  c3 

q(l,  3)  =  c2:  q(2,  4)  =  c2:  q(3,  1)  =  c2:  q(4,  2)  =  c2 
el  =  Application. WorksheetFunction.MMult(Phi,  Sigma) 
e2  =  Application.WorksheetFunction.Transpose(Phi) 
e3  =  Application.WorksheetFunction.MMult(el,  e2) 

Mu  =  Application.WorksheetFunction.MMult(Phi,  Mu) 

Call  matadd(e3,  q.  Sigma) 

’Defining  Mu_W  (long-term  average  velocity)  in  terms 
'of  horizontal  east-west  velocity  (HorV) 

'and  in  terms  of  vertical  north-south  velocity  (VerV) 

Dim  HorV#,  VerV# 

'To  change  the  long-term  average  velocity  from  zero  to 

'an  input  value  in  cells  s  1 3  and  si 4,  comment  the  next  two  lines  and  uncomment 
'the  lines  defining  HorV  and  VerV  as  called  values. 

HorV  =  0 
VerV  =  0 

'HorV  =  Worksheets("master").Range("sl3").Value 
'VerV  =  Worksheets("master").Range("sl4").Value 
'Adding  Mu_W  for  the  long-tenn  average  velocity  of  the  target 
Mu(3,  1)  =  Mu(3,  1)  +  HorV  *  (1  -  b2) 

Mu(4,  1)  =  Mu(4,  1)  +  VerV  *  (1  -  b2) 

End  Sub 


7.  Orientation  of  Compass 

Public  Function  Atn2(x,  y) 

’-pi<=atn2<pi 

'counterclockwise  from  East  in  radians 
’Atn2(y,x)  is  clockwise  from  North  in  radians 
Const  pi2  =  1.570796326794 
If  x  >  0  Then 
Atn2  =  Atn(y  /  x) 

Elself  x  <  0  Then 
If  y  >  0  Then 

Atn2  =  Atn(y  /  x)  +  pi2  +  pi2 
Else 

Atn2  =  Atn(y  /  x)  -  pi2  -  pi2 
End  If 

Elself  y  >  0  Then 
Atn2  =  pi2 
Else 

Atn2  =  -pi2 
End  If 
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End  Function 


8.  Course  and  Bearing  Functions 

Function  bearing(latl,  lonl,  lat2,  lon2) 

'course  in  degrees  from  N  clockwise,  inputs  in  degrees,  Wlon<0 
Dim  conv#,  xl#,  x2#,  x3#,  x4# 
conv  =  Atn(l)  /  45 

xl  =  latl  *  conv:  x2  =  lonl  *  conv:  x3  =  lat2  *  conv:  x4  =  lon2  *  conv 
bearing  =  course(xl,  x2,  x3,  x4)  /  conv 
End  Function 

Function  course(latl,  lonl,  lat2,  lon2) 

'course  in  radians  from  N  clockwise  to  get  from  1  to  2 
Dim  cd#,  dlon#,  cl#,  c2#,  si#,  s2# 
cl  =  Cos(latl) 

Ifcl  <=  0  Then  'north  or  south  pole 
If  latl  >  0  Then 
course  =  4  *  Atn(l) 

Else 

course  =  0 
End  If 
Else 

dlon  =  lon2  -  lonl 

c2  =  Cos(lat2):  si  =  Sin(latl):  s2  =  Sin(lat2) 

cd  =  si  *  s2  +  cl  *  c2  *  Cos(dlon) 

course  =  Atn2((s2  -  si  *  cd)  /  cl,  c2  *  Sin(dlon)) 

If  course  <  0  Then 

course  =  course  +  8  *  Atn(l) 

End  If 
End  If 

End  Function 


9.  Matrix  Algebra 

Sub  matadd(A,  B,  c) 

Dim  i&,  j& 

For  i  =  1  To  UBound(A,  1) 

For  j  =  1  To  UBound(A,  2) 
c(i,j)  =  A(i,j)  +  B(i,j) 

Next  j 
Next  i 
End  Sub 

Sub  matdif(A,  B,  c) 

Dim  i&,  j& 

For  i  =  1  To  UBound(A,  1) 
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For  j  =  1  To  UBound(A,  2) 
c(i,j)  =  A(i,j)  -B(i,j) 

Next  j 
Next  i 
End  Sub 

Sub  matscal(A,  B,  c) 

Dim  i&,  j& 

For  i  =  1  To  UBound(A,  1) 

For  j  =  1  To  UBound(A,  2) 
c(i,j)  =  A(i,j)  *  B 
Next  j 
Next  i 
End  Sub 

B.  UPDATING  AOU  TO  FUTURE  TIME  T 

Private  Sub  UpdateTo_Click() 

With  Worksheets("targetl") 

Sigma  =  .Range("al:d4").Value 
Mu  =  .Range("a6:a9").  Value 
told=  .Range("CurTime")(l,  l).Value 
tnew  =  .Range("CurTime")(2,  l).Value 
If  told  <  tnew  Then 

Call  KMove(Mu,  Sigma,  tnew  -  told) 
.Range("al:d4").Value  =  Sigma 
.Range("a6:a9"). Value  =  Mu 
.Range("CurTime")(l,  1).  Value  =  tnew 
.Range("CurTime")(2,  l).Value  =  "" 

End  If 
End  With 
End  Sub 


47 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


48 


BIBLIOGRAPHY 


Bell,  Bradley  M.,  and  Fredrick  W.  Cathey.  “The  Iterated  Kalman  Filter  Update  as  a 
Gauss-Newton  Method,”  IEEE  Transactions  on  Automatic  Control  38,  no.  2 
(February  1993):  294-297. 

Brown,  Robert  G.,  and  Patrick  Y.  C.  Hwang.  Introduction  to  Random  Signals  and 

Applied  Kalman  Filtering:  with  Matlab  exercises  and  solutions.  (3ld  Edition)  New 
York:  Wiley  &  Sons,  1997. 

Clark,  Vem.  Anti-Submarine  Warfare  Concept  of  Operations  ( Draft  26  April  2004). 

Czech,  Carl.  U.S.  Navy  Publication.  NPRDC-TR-98-2  The  Interactive  Multisensor 

Analysis  Training  (IMAT)  System:  an  Evaluation  of  Airborne  Acoustic  Mission 
Course.  San  Diego,  CA:  GPO,  1998. 

Fukumori,  Ichiro.  “A  Partitioned  Kalman  Filter  and  Smoothers,”  Monthly  Weather 
Review  130  (May  2002):  1370. 

Kerr,  Thomas  H..  “Streamlining  Measurement  Iteration  for  EKF  Target  Tracking,”  IEEE 
Transactions  on  Aerospace  and  Electronic  Systems  27,  no.  2  (March  1991):  408- 
421. 

Maybank,  S.  J..  “Bearings-Only  Tracking  in  the  Plane,”  SIAM  Journal  on  Applied 
Mathematics  58,  no.  3  (June  1998):  975-998. 

Olson,  Mark  A..  “Simulation  of  a  Multitarget,  Multisensor,  Track-Splitting  Tracker  for 
Maritime  Surveillance,”  Naval  Postgraduate  School  Master’s  Thesis,  (September 
1999). 

Robbins,  Douglas  Leigh.  “Decision-Making  Process  of  an  Antisubmarine  Warfare 
Commander”  Naval  Postgraduate  School  Master 's  Thesis,  (September  1986). 

Sen,  Adnan,  Hakan  A.  Cirpan,  and  Erdal  Panayirci.  “Joint  Channel  Tracking  and  Symbol 
Detection  for  OFDM  Systems  with  Kalman  Filtering,”  AEU  International  Journal 
of  Electronics  and  Communications  57,  no.  5  (2003):  317-327. 

U.S.  Navy  Publication.  NWP  3-21.51.3  Surface  Ship  Passive  Localization  and  TMA. 
Washington,  D.C.:  U.S.  GPO. 

U.S.  Navy  Publication.  NTTP  3-21.23  Submarine  Tracking  Manual.  Washington,  D.C.: 
U.S.  GPO. 

49 


U.S.  Navy  Publication.  NWP  3-10.3  Inshore  Antisubmarine  Warfare.  Washington,  D.C.: 
U.S.  GPO. 

U.S.  Navy  Publication.  NWP  3-21  Navy  ASW.  Washington,  D.C.:  U.S.  GPO. 

U.S.  Navy  Publication.  NWP  3-21.51.1  Target  Motion  Analysis.  Washington,  D.C.:  U.S. 
GPO. 

Wagner,  Daniel  H..  “Naval  Tactical  Decision  Aids,”  Military  Operations  Research 
Lecture  Notes,  (September  1989). 

Washburn,  Alan.  A  Short  Introduction  to  Kalman  Filters.  NPS:  2004. 

Woodward,  John.  One  Hundred  Days:  the  Memoirs  of  the  Falklands  Battle  Group 
Commander.  Annapolis:  Naval  Institute  Press,  1992. 


50 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Technical  Infonnation  Center 
Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 
Naval  Postgraduate  School 
Monterey,  California 

3.  Professor  Washburn 
Naval  Postgraduate  School 
Monterey,  California 

4.  Vice  Admiral  Bacon,  USN  (Ret.) 

Naval  Postgraduate  School 
Monterey,  California 

5.  Captain  Ferguson,  USN 
Carrier  Group  One 
San  Diego,  California 

6.  N-75,  Director  Expeditionary  Warfare  Division 
Chief  of  Naval  Operations 

Washington,  D.C. 

7.  N-77,  Director  Submarine  Warfare  Division 

Chief  of  Naval  Operations 
Washington,  D.C. 

8.  Fleet  ASW  Command 
San  Diego,  California 


51 


9.  Commander  Submarine  Development  Squadron  Twelve 
Naval  Submarine  Base 

Groton,  Connecticut 

10.  Commander  Third  Fleet 

ATTN  Naval  Postgraduate  School  Representative 
San  Diego,  California 

1 1 .  Commander  Submarine  Group  Seven 
Yokosuka,  Japan 

12.  Commander  Submarine  Group  Eight 
Naples,  Italy 

13.  Office  of  Naval  Research 

Ocean,  Atmosphere  and  Space  Department 
Arlington,  Virginia 

14.  Commander  U.S.  Pacific  Fleet 

ATTN  N00/ASW  Captain  Brandhuber,  USN  (Ret.) 
Pearl  Harbor,  Hawaii 


52 


